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Abstract 

The relationships between thermal conductivity and other petrophysical properties 
have been analysed for a borehole drilled in a Tertiary Flysch sequence. We establish 
equations that permit us to predict rock thermal conductivity from logging data. A 
regression analysis of thermal conductivity, bulk density, and sonic velocity yields 
thermal conductivity with an average accuracy of better than 0.2W(mK)-^ As 
a second step, logging data is used to compute a lithological depth profile, which 
in turn is used to calculate a thermal conductivity profile. From a comparison of 
the conductivity-depth profile and the laboratory data it can be concluded that 
thermal conductivity can be computed with an accuracy of less than 0.3 W(mK)"^ 
from conventional wireline data. The comparison of two different models shows that 
this approach can be practical even if old and incomplete logging data is used. The 
results can be used to infer thermal conductivity for boreholes without appropriate 
core data that are drilled in a similar geological setting. 
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1 Introduction 

The German Molasse Basin is currently analysed with geothermal methods in 
order to determine groundwater flow rates in the deep subsurface ([ll). This 
requires to determine accurately and separate from each other the effects 
of the different heat transport processes: Steady-state conductive, transient 
conductive, and advective heat flow. Computing the steady-state conductive 
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temperature field requires detailed thermal conductivity data. Two require- 
ments should be met: First, lateral variations in thermal conductivity have to 
be known accurately for regional studies. Second, a detailed analysis of a par- 
ticular borehole requires a high resolution thermal conductivity profile. Core 
material in that area is either confined to a few boreholes or available only 
for a particular layer. This violates both requirements. Logging data, on the 
other hand, offer better spatial coverage and depth resolution. In a situation 
with an incomplete coring sequence this helps to prevent biased sampling in 
lithologies not representative of the full sequence. It is therefore highly de- 
sirable to derive thermal properties from logging data to improve the sparse 
geothermal database. 

Methods to compute thermal conductivity from wireline logs can be classified 
into two main categories (jij). The first approach relates one or more logging 
measurements or some derived property directly to thermal properties via 
empirical relationships. This method has been used to compute thermal con- 
ductivity in several studies (jsj; 0; S @) and a review is given in (0). In the 
second approach, the major mineral or rock components are identified and 
the volumetric fractions of these components are derived from regular wire- 
line data. This composition together with component thermal conductivity 
values is then used to compute the effective thermal conductivity assuming an 
appropriate mixing law (0; [ol: fiol). This approach is more flexible than the first 
one. For instance, it does not require the same suite of logs in each borehole 
as long as they are suitable to compute component volumes. 

It has been pointed out (0; Q) that results of the current methods are confined 
to a geographic region or geological setting. In the following we would like to 
assess the two different methods discussed above for the Molasse Basin. As 
the area was explored for hydrocarbons in the past, logging data often are old 
and comprise only sonic or natural gamma logs. Although it would be more 
desirable to use only the second approach, it is necessary to evaluate direct 
methods as well to make full use of the available logging data. We test direct 
methods on a set of core data, that was analysed in the laboratory and the 
compositional approach using the same set of laboratory data together with 
wireline data obtained in the borehole. 



2 Core data 

The data set for laboratory analysis is obtained from a borehole drilled for 
hot water into the Upper Marine Molasse formation. This formation consists 
of shaly sandstones and marls. The borehole was cored between 570-810 m 
depth with nearly 100 % recovery. The sequence consists of a succession of 
shaly sandstones and marl beds. Porosity ranges from 10-30 %. Split cores 
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were available for measurements in sections of 1 m length. In general, we anal- 
ysed every fifth section, but in some depth intervals we studied all cores. The 
cores had broken into several pieces during coring and splitting. We selected 
two to five samples from each core section for further analysis. The samples 
were dried at 60°C to prevent cracking or alteration of clay minerals. We 
measured thermal conductivity, sonic velocity, and bulk density on the dried 
samples. Then the samples were saturated under vacuum and the same mea- 
surements were repeated for the saturated samples. Because the samples were 
not well consolidated a number of samples was damaged or destroyed during 
saturation, resulting in a lower number of saturated measurements. Figure [T] 
shows histograms of the results. 



Thermal conductivity was measured using the optical scanning technique 
which yields a continuous profile of thermal conductivity along the core axis of 
the sample (jlll : Il2l : Il3l ). We also measured thermal conductivity perpendicu- 
lar to the core axis to obtain anisotropy. Anisotropy ratios were generally less 
than 3-4% so that we assume the rock to be isotropic. Values for each sample 
were computed as averages of the scanning line. We measured bulk density and 
sonic velocity with a multi-sensor core- logger (e. g. The measurement is 

performed perpendicular to the core axis. Both bulk density and sonic velocity 
records need special attention before they can be used in the analysis. Bulk 
density is measured by gamma-ray absorption mainly due to Compton scatter- 
ing (jl5l ). It is calibrated routinely using an aluminium standard. Since electron 
density is measured rather than real density, better results are obtained when 
using a calibration standard with a similar electron density (jl5l ). Therefore 
we measured bulk density and porosity independently on a number samples 
using either Archimedes' method or gas-state and solid-state pyknometers. 
The obtained bulk density is then plotted versus the absorption coefficient fi 
measured with the core-logger for dry and saturated samples. Figure [2] shows 
that a linear relationship exists which can be used for calibration. For dry and 
saturated measurements the following equations are obtained: 



Pfe_s = 161.5yUs (saturated) (1) 
Pb4 = l62.1fLd (dry) 

The slope is slightly lower for saturated measurements. A detailed analysis 
shows that saturated gamma density should be corrected using a separate 



porosity measurement (jl6l ). However, considering the data scatter and the 
limited range of the saturated measurements this can be included in the single 
calibration factor. Sonic velocity was measured on the samples under ambient 
pressure and we found a large difference in velocity obtained from wireline logs 
and cores (Figure [3l bottom panel). It is well known that changes in pressure 



have a profound influence on elastic properties (e. g. (1171 : llSl)). An empirical 



relationship with linear and exponential terms is usually used to correct this 
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effect. In order to quantify tlie correction, sonic velocity was measured on a 



sample under uniaxial pressure as described in (llTl ). We fitted two exponential 
equations to this data (Figure [3]), both with an exponential term but one 
without a linear term . The following coefficients were obtained, with pressure 
P in MPa and velocity v in ms"^: 

v{p) = 2360 + 787e°-356-P +35. 7P (2) 
^;(p) = 2760 + 1180 e°-24=^^ (3) 

The rms-errors are 19 ms"^ (eq. [2]) and 18 ms"^ (eq. [3]). Both equations fit the 
pressure-velocity data equally well and differ only in the extrapolated range 
from 450 m to 820 m. This is, however, the depth range of the samples and 
equation [3] yields better results when we compare corrected laboratory data 
to logging data. Equation [2] slightly over-corrects the laboratory data. This 
is probably caused by the fact that the measurements only cover a range of 
up to 10 MPa. In this pressure range the exponential term is more significant, 
and the linear term is not well defined. Therefore, we used equation [3] to 
correct the measurements (Figure [31 lower part). The remaining scatter after 
correction may be due to the use of only one set of fit parameters. A detailed 
analysis should take into account the differing composition of the samples. 
Because of these limitations, we did not use sonic velocity in our analysis 
of laboratory and logging properties. We did use it, however, for correlation 
analysis of properties measured in the laboratory. Following this preliminary 
work of preprocessing laboratory data the next step is to examine correlations 
of properties measured in the laboratory. We discuss the results of this work 
in the next section. 



3 Laboratory results 



Figure [T] shows a summary of the measurements processed as described in the 
previous section. Dry properties show a broader distribution of values, in the 
case of thermal conductivity it is even bimodal. This follows from the fact 
that the ratio of fluid and matrix properties is larger for the rock/air-system 
than for the rock/water-system. Hence, also the range of the effective values 
is larger for dry properties. 

Our goal in measuring the various properties is to predict thermal conductivity 
from other petrophysical properties which can be measured in-situ. For this 
purpose we analysed linear correlations between thermal conductivity in dry 
and saturated condition and and other measurements in their corresponding 
saturation state. Figure |4] shows cross-plots of these data. Porosity was com- 
puted from dry and saturated bulk density. We did not use sonic measurements 
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to compute porosity as velocity-porosity models work best when velocity is 
measured under confining pressure such that the rock is at its terminal velocity 



(Il9l: 1201) ■ We grouped the samples according to their lithology: The first group 
contains shaly sandstone samples ("sandy samples"), the second group holds 
marlstones ("carbonaceous samples"). This division is somewhat preliminary 
as it is based on geological core descriptions and not on an analysis of the 
mineral assemblages. 

A regression analysis of thermal conductivity versus sonic velocity, bulk den- 
sity and porosity, respectively, was performed. For each combination a linear 
equation of the form 

y = aix + Oq (4) 



was fitted to the data. The regression analysis uses a total-least-squares solu- 
tion to account for measurement errors both in the dependent and independent 
variables (12 ll) . The Jackknife- method yields the variances of the parameter 
estimates fl22l ). We also coraputed correlation coefficient R and rms-error to 



assess the quality of the fit (1231 ) 



Table [T] summarises the results. If the rms-error is interpreted as the predictive 
error of the computed relationships it can be deduced that it is possible to 
calculate thermal conductivity from density and or sonic velocity with an 
average accuracy of about 0.1-0.2 W(mK)"^, given an appropriate data set for 
calibration. Correlation coefficients are largest for dry properties. This is due to 
the larger contrast in rock and pore volume properties for dry samples. Thus, 
porosity variations cause larger variations in the effective properties of the two- 
component system and therefore stronger correlations for dry samples. On the 
other hand, the low correlation coefficient between thermal conductivity and 
porosity is due to the complex relationship between these properties. Whereas 
a linear equation will be sufficient to describe the relationship between thermal 
conductivity and other petrophysical properties, it is inadequate for porosity. 
This is addressed in more detail in section |H Grouping samples according to 
their lithology has a larger effect on the quality of the fit for saturated than for 
dry samples. Figure H] shows in the crossplot of thermal conductivity versus 
sonic velocity or bulk density that the two lithology groups have the same 
general trend for dry samples but are separated for saturated measurements. 
This again is an effect of the strong influence of porosity for dry samples which 
masks variations due to lithology. From our analysis it thus appears that dry 
measurements are most sensitive to porosity changes whereas saturated ones 
reflect both variations in lithology as well as porosity. 

The question arises if a correlation using more than one petrophysical property 
improvea the prediction of thermal conductivity. Table [2] shows the results of 
a multiple linear regression for thermal conductivity. The equation has the 
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form: 



A = ao + aiVp + + a3</> (5) 



There is only slight improvement over the simple regressions. In practise this 
improvement might be outweighed by the additional effort for performing the 
measurements. If the samples are considered to be composed of the three 
components sand, shale, and carbonate, the quality of the predicted thermal 
conductivity is directly related to the quality of prediction of these three com- 
ponents from the measured properties. However, quartz and calcite differ in 
density and velocity by only 3 and 13 %, respectively. Shales have a large 
variability in their physical properties, but they usually differ considerably 



from quartz and calcite (jl5l : |2J). Therefore, any of the measurements will be 
more sensitive to variations of shale content and porosity than to changes in 
carbonate and sand content. Thus, a combination of sonic velocity and bulk 
density does not provide significantly more information than each of them 
alone. This fact can be assessed in a cross-plot of these two properties (Fig- 
ure HJ lower left panel). Sandy and carbonaceous samples essentially plot on 
top of each other. Additional measurements, such as natural gamma radiation 
for instance, are required to characterise our samples better. The contrast in 
thermal conductivity between quartz and calcite is about 60 %, sufficiently 
large to separate them in a cross-plot (Figure HI upper left). This offers oppor- 
tunities for characterising lithology, as conductivity can be determined rapidly 
and continuously along a core. 



4 Choosing an appropriate mixing law 



A mixing law for calculation of the effective thermal conductivity of a compos- 
ite medium according to the content and thermal conductivity of its compo- 
nents is required both in the analysis of thermal conductivity measured in the 
laboratory on dry and saturated samples as well as in the prediction of ther- 
mal conductivity from logging data. Several models have been proposed and it 
is instructive to review the relationships most commonly used in geothermics 
(0). Effective thermal conductivity A of a layered medium with thermal con- 
ductivities Ai and A2 depends on the direction of the temperature gradient. If 
heat flow is parallel to the layering the effective conductivity is equal to the 
arithmetic mean (Aa) layer thermal conductivities, weighted by their volume 
fractions. For perpendicular heat flow it corresponds to the harmonic mean 



(A,) my- 



Xa = Vi\i + V2\2 (6) 



6 



While these averages are useful to estimate the average thermal conductiv- 
ity of a vertical rock sequence they are inappropriate for estimating effective 
sample thermal conductivities. Narrower bounds can be derived by assuming 
a geometry where the solid consists of spheres dispersed in the pore fluid or 



were the fluid is confined in spherical inclusions in the rock matrix (l26l : |2 



This configurations yield the lower and upper Hashin-Shtrikman (HS , HS^ 
bounds, respectively: 



Ahs- - Ap + 1 ^ ^ - (8) 

Ahs+ = Am H ^ — ^ (9) 

Here denotes porosity, and the subscripts p, m denote pore and rock matrix 
properties, respectively. These bounds are of theoretical importance because 
effective thermal conductivities of rock samples should generally fall in be- 
tween these bounds. However, in many cases they are to far apart to be of 
practical use. In this situation an estimate can onl y b e obtained from empir- 



ical relationships such as the geometric mean (1281 ) . that is often used in 
geothermal studies: 

A, = XtXlr^'^ (10) 



However, other researcher preferred to use the average of the upper and lower 
HS bounds (1271 ) , or the square root average (1291 : (2) : 



Xs = (f)JXp + (1 



'111 



The self consistent approach (1301 : l3ll ) is popular for elastic properties but is 
not widely used in geothermal research. For porosities typical for rocks it gives 
results similar to the square-root average. The thermal conductivity A;, for a 
two component medium is given by the equation: 



2/3 + Ap/(3Ab) 2/3 + Xm/iSXt) 



:i2) 



The particular choice of a model becomes important when the contrast in 
thermal conductivity of the constituents increases. Figure [5] illustrates this by 
showing bounds and estimates for the saturated and the dry case. Thermal 
conductivity is assumed to be 5W(mK)"^ for the matrix and 0.6W(mK)"^ 
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and 0.026 W(mK)"^ for water and air, respectively. In this case the maximum 
difference between lower and upper HS bounds are 3.7 and l.OW(mK)"^ for 
dry and saturated samples. For a two-phase mineral assemblage of crystalline 
quartz (7.7W(mK)-i) and orthoclase (2.3W(mK)-i) the maximum difference 
is as low as 0.4W(mK)"^. This shows that the choice of a correct mixing law 
for a mineral assemblage is somewhat arbitrary while it is essential when con- 
sidering dry samples. Although a dry sample will rarely occur in nature, these 
considerations are important when laboratory measurements on dry samples 
are used to predict in-situ saturated thermal conductivity. 

Another important conclusion is, that in the case of a mineral assemblage the 
geometric mixing law closely follows the lower Hashin-Shtrikman bound, corre- 
sponding to a rock model consisting of grains suspended in a fluid. The square 
root law, on the other hand, is very close to the upper Hashin-Shtrikman bound 
and could be interpreted as a well lithified rock with spherical pores. Thus, 
as each particular empirical mixing laws corresponds to a particular rock- 
structure, one single model cannot be adequate for all different rock struc- 
tures. 

Additional parameters can be introduced in order to incorporate rock structure 
into a mixing law. Several models assume spheroidal pores where a is the 



aspect ratio of the spheroids (1321 : l33l : l34l : l35l ) . In an application of the model 



by Zimmerman (1331 ) aspect ratios as low as 0.1 were found for basalts (1361 ). 
much less than the actual aspect of the pores. It was regarded as a value 
representing the aspect ratio of the grain contact rather than that of the 
pores. 

From this discussion it appears that mixing laws of varying complexity can be 
used in different situations: Thermal conductivity of a mineral assemblage or 
a water-saturated rock can be computed with sufficient accuracy using a geo- 
metric mixing law. For measurement on dry samples the pore structure needs 



to be taken into account. The model proposed in (1331 ) was therefore applied to 
our laboratory data set in order to determine the aspect ratio a. The model 
assumes a homogeneous mixture of randomly distributed spheroids. For a rock 
with oblate spheroidal pores one obtains the effective thermal conductivity A: 

A (1 -0)(1 -r) -fr/5(/) 



Am (l-0)(l-r)+/3</. 
The parameters /3, M, and Q are defined by: 



(13) 
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M = 

e= 



3 \2 + {r-l)M l + (r-l)(l-M) 
29 - sin 29 
2tan9 8111^9 



arccos a 



(16) 
(17) 



Again, a is the aspect ratio of the spheroidal inclusions. The model described 
by equations [13] - [T7| consists of five parameters: effective, matrix, and pore 
fluid thermal conductivity, porosity, and aspect ratio. Our measurements of 
saturated and dry thermal conductivity yield two equations of this type which 
differ only in the effective thermal conductivity A and the pore fluid thermal 
conductivity A/. This yields seven parameters in total. If we use the two ther- 
mal conductivity measurements, known values of water and air thermal con- 
ductivity, and an independent porosity measurement, we obtain two equations 
with two unknowns: Matrix thermal conductivity and pore aspect ratio. We 
solved the system for our data by a nonlinear iterative algorithm. Figure M 
shows a crossplot of the computed values of matrix thermal conductivity A^ 
and aspect ratio a. The median values of a are 0.011 and 0.016 for sandy and 
carbonaceous samples, respectively. This is about one magnitude less than 
the values reported in (f36l ). However, those values were measured on igneous 
rocks and crack aspect ratios based on sonic or mechanical measurements gen- 
erally range from 10-2-10"^ (i33;l3i). Median matrix thermal conductivities are 
5.14 W(mK)"^ and 4.09 W(mK)"^ for sandy and carbonaceous samples. 



5 Thermal conductivity predicted from wireline data 



Wireline logs of natural gamma radiation (GR), neutron porosity (NPHI), and 
bulk density (RHOB) were available for analysis in the depth range 570-800 m. 
In the depth range from 600-800 m caliper (CALI) and temperature (TEMP) 
had been logged during an aquifer test. The logs are shown together with the 
core lithology in figure [71 The temperature log is strongly disturbed by water 
flowing from the formation with temperatures varying from 35-37°C. Unfor- 
tunately the disturbance renders a quantitative interpretation of temperatures 
with regard to the conductive regime impossible. The hole was drilled with 
8.5" diameter and the caliper generally reads below 9". Some larger diameter 
sections (630-680 m) are consistent with lithological changes. The core lithol- 
ogy showed no signs of hydrocarbons and the drilled sequence is not known 
as a hydrocarbon reservoir. Thus, log data quality can be regarded sufficient 
for quantitative interpretation. 

We applied several editing steps before we analysed the log ging data. For a 



general discussion of these steps see for example (l39l : l4d : Il5l ) . Logging curves 
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from different tool runs were corrected to common depth points, bad data were 
eliminated, and environmental corrections applied. Core depths are shifted to 
match logging depth. Core data are smoothed before they are compared to 
log curves since these have a lower depth resolution. For this purpose we used 
an inverse distance algorithm that employs an Gaussian weighting function: 



N 



Wk = Y.^~^^' (18) 
1=1 

— ^Wie"(^)' (19) 



i=l 
1 ^ 

Vk 



Here are data values, Zi is the depth of the data value, and Wk is the 
normalisation constant. We took the weighting distance to be 0.5 m. 

In general the response W of a wireline tool j is determined by the volume 
fractions Vi of the rock components i and their theoretical log response T/. 
Assuming a linear relationship this results in 

w = Y.v{ri (20) 



with the constraint that I]i = 1 (14 ll : |42| ) . This equation is correct for prop 



erties like bulk density or neutron porosity. For acoustic properties, on the 
other hand, it is not generally valid but results in empirical relationships. 
If for instance the slowness is considered as the sonic measurement, equa- 



tion [20] implies Wyllie's classical travel-time average (1431 ). If the number of 
constituents equals the number of tool responses equation [20] has one solution. 
If the number of equations is larger than the number of components, the sys- 
tem is overdetermined and equation [2D] can be solved in a least-squares sense. 
Thus, lithological composition can be computed from a sufficient number of 
logs. Based on this composition and the known values of thermal conductivity 
of the components, an effective thermal conductivity can be computed. The 
geometric average law is used to compute the effective thermal conductivity 
of the mineral mixture. This value will be used as the matrix thermal con- 
ductivity. The effective conductivity of dry and water saturated rocks is then 
calculated from equation [T3] using a median aspect ratio of 0.012. As we want 
to compare the results to our laboratory data, no temperature correction is 
necessary at this stage. 

As discussed previously for laboratory data, we explore different models with 
varying simplicity to evaluate how well thermal conductivity can be described 
by logging data under different circumstances. Two models are of particular 
interest: (1) One model consists of a mixture of sand, shale, and carbonate and 
uses the information of all logs available for the borehole; (2) Another model 
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consists only of the two components sand and shale, ignores the carbonate 
fraction, and uses only wirehne logs of slowness and natural gamma radiation, 
DT and GR, respectively. This model is required for older wells where only 
these two logs are available. Although these wells have not been logged with 
modern tools they comprise a large fraction of the available dataset. 

The full model uses logs of gamma density (RHOB), natural gamma radiation 
(GR), and slowness (DT) as input logs. We employed a commercial software 



package specialised in deriving rock composition from wireline logs (144| ). The 
three components sand, shale, and carbonate were parameterised in the inver- 
sion using response values of the minerals quartz, glauconite, and calcite. The 
choice of the shale mineral glauconite, an iron-rich variety of illite, is based 
on the geological description of the cores and general information about the 



geology of the Upper Freshwater Molasse (1451 ). It might not be the only shale 
mineral present, but the values for illite, another abundant shale mineral, are 
very similar to the ones we used. Values for the logging properties and min- 
eral thermal conductivities are summarised in Table [31 Thermal conductivity 
values for quartz and calcite are 7.69 W(mK)"^ and 3.59 W(mK)"^, respec- 
tively (j46l ). The choice for the shale fraction is more difficult. Glauconite was 
measured with a value of 1.6W(mK)"^ (j27l ). We found that an optimal fit 
can be obtained with a higher value of 2.2W(mK)"^. This discrepancy could 
be easily explained by the variability of properties for shale minerals, but 
could also indicate that small amounts of other minerals are present which 
are not accounted for. There is generally a good agreement between computed 
and measured thermal conductivity (Figure [8]). The rms error of the recon- 
struction is 0.27 W(mK)"^ and 0.28 W(mK)"^ for saturated and dry thermal 
conductivity. Although the rms-error is slightly larger for dry properties, they 
can be reproduced much better than saturated properties because the large 
conductivity contrast of air/rock matrix enhances the variations in thermal 
conductivity. 

The second predictive model employs a simple mineralogy consisting only 
of sand and shale. Additionally, an empirical relationship is used between 
sonic velocity Wp, shale content Vshaie, and porosity (j), established for shaly 



sandstones and a pressure range from 10-40 MPa (1471 ). At a pressure of 20 MPa 



(about 900 m) the equation for the sonic velocity is given by 

Vp[km s-i] = 5.49 - 6.940 - 2.17Vshaie. (21) 



The empirical parameters in this equation are slowly varying functions of 
pressure. For the gamma-ray log we use the linear log response equation [20] 
assuming no radiation for the fiuid: 

GR = GTisandVsand + GRs haleVs hale (22) 
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The parameters used for the model are given in Table [31 Using the mean 
density of our laboratory measurements of 2330 kgm"^, depths were converted 
to lithostatic pressure which is needed for equation [2TJ The resulting lithology 
profile is shown in figure M 

Again we test this model for its ability to predict thermal conductivities. For 
this purpose we sampled the compositional log at the depths of the core mea- 
surements. Computed composition and measured values of thermal conduc- 
tivity were then used in an inversion to find optimal mineral thermal conduc- 
tivities of the constituents (Table [3]). A synthetic log of thermal conductivity 
is then computed and compared to the measurements. The rms-misfit of this 
two-component model is 0.28 W(mK)"^, essentially the same as the misfit of 
the three-component model. 

When comparing the lithology logs for the two models it is apparent that 
the volume fraction of the shale component has not changed very much. The 
sand fraction of the two-component model also includes the missing carbonate 
fraction. This is also reflected in the lower mineral thermal conductivity of the 
sand/shale model. As a consequence the model will be only successful when 
the ratio of sand and carbonates does not change too much. This can be seen 
in figure [9] at depths around 630 m, where high thermal conductivities cannot 
be accurately reproduced by the model. 

Both models display a better fit in the upper part of the profile than in the 
lower part below 700 m. We believe that this is due to an overestimation of 
porosity for the cores we analysed. Average porosity in the depth range 750- 
800 m is 20% for core data but only 17% for log data. The difference can be 



attributed to the release of overburden pressure during coring (|48| ). However, 
at this point no measurements under confining pressure could be performed to 
verify this effect. Depending on the geological setting and method to compute 
an overburden correction, correction factors of 0.85-0.95 have been reported 



(l49l) . It is to be expected that our samples respond strongly to the pressure 
relief due to their poor consolidation. In this situation log derived porosities 
might be more accurate and reliable than core derived porosities adding to 
the usefulness of log derived thermal conductivity. 

So far, conductivities were computed only for room temperature. For model- 
ing of temperatures in-situ it is necessary to correct the effect of temperature 
on thermal conductivity. Temperature dependent thermal conductivity was 
measured for 23 samples of the Molasse Basin, one of them from the borehole 
analysed here (jH). An empirical function based on (jsol ) of the form 
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was used fit the data (15 ll ). Here Aq is the thermal conductivity at room 
temperature (25°C). Coefficients are a = 0.99, b = 3.4 ■ 10~^ K"^, and c = 
3.9 • 10"'^ W(mK)"^. Using the aquifer temperature of 36°C at 800m depth 
and a range of 2.5-3.5 W(mK)"^ for saturated thermal conductivity the nec- 
essary correction amounts to 5-7%. 



6 Summary and Conclusion 

Empirical relationships between thermal conductivity and other petrophysical 
properties depend on local conditions, in particular the type of diagenesis for 
the rocks. We analysed core and log data from one borehole in the Molasse 
Basin in order to establish a set of such equations. Empirical relationships 
were derived from sonic velocity or bulk density laboratory data that allow 
to predict thermal conductivity to an accuracy of about 0.2W(mK)"^, on 
average. Predicting thermal conductivity from logging data has a larger error 
of about 0.3 W(mK)"^. On the one hand, this may be expected in view of the 
different spatial resolutions and the problems encountered in the matching 
of core and log depths. On the other hand, this is outweighed by the larger 
number of sampling points. An important restriction, as with other studies of 
this type, is that the results are restricted to the particular conditions in a 
specific basin. 

An ideal combination of wireline logs for an optimum determination of lithol- 
ogy would comprise the full suite of nuclear measurements. In case of the Mo- 
lasse Basin, we demonstrate that thermal conductivity can be derived even 
though the number of available logs is less than ideal. This is an important 
result, as for many old wells - drilled before the advent of modern logging 
tools - often only a natural gamma log is available. The fact that thermal 
conductivities can be reconstructed from these data alone makes these older 
wells attractive for such an analysis. 

The discussion of laboratory data and the two compositional models shows 
that this log combination is most useful in a setting dominated by sand and 
shale fractions. In contrast, a combination of sandstone and carbonate cannot 
be well characterised. Fortunately though, these two components differ by a 
large contrast in thermal conductivity. This measurement can be performed 
rapidly and yields a continuous proffie along the core. Also, thermal conduc- 
tivity is closely linked to other petrophysical properties, such as permeability 
(jiil). This can be particularly useful, for instance, in reservoir studies where 
permeability or porosity are important properties. 

One aspect which we did not address in our analysis is the variation of thermal 
conductivity with in situ pressure. Sonic velocity varies strongly with pressure 
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and so does thermal conductivity (1461: |52|) . Attempts have been made to model 
this behaviour with respect to the variation of elastic properties with pressure 
( I53I : l33l ). Because the Tertiary samples examined in our work are poorly con- 
solidated, a considerable pressure effect may be expected. This requires more 
detailed study in future. 
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Table 1 

Results of single regression analysis of thermal conductivity A of core samples based 
on sonic velocity Vp, bulk density pb, and porosity cp. The top panel gives results for 
dry, the lower part for saturated samples. For each regression the fit parameters ai 
and oq (Equation m , with errors, the correlation coefficient R, and the rms-error of 
the fit are shown. Fits were computed both for all samples and grouped by lithology. 
For ease of display units of Vp and pi, are kms"^ and gcm"'^, respectively. 
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Table 2 



Results of the multiple regression of thermal conductivity A versus sonic velocity 
Vp, bulk density pb, and porosity (p using equation [5j For each regression coefficients 
of the fit, coefficient error, correlation coefficient R, and the rms-error of the fit 
are given. Fits were computed both for all samples (lower part) and grouped by 
lithology (upper part). For ease of display units of Vp and pt are kms"^ and gcm"'^. 
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Table 3 

Summary of response parameters and component thermal conductivities used in the 
log data analysis. Slowness DT, natural gamma radiation GR, and neutron porosity 
NPHI are used in the complex model. For the simple model NPHI is not used and 
DT coefficients are determined from equation 1211 
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Fig. 1. Histograms of sonic velocity Vp, bulk density ph, and thermal conductivity A 
measured in the laboratory on core samples. Properties measured are from top to 
bottom: sonic velocity, bulk density, and thermal conductivity. Values are shown for 
saturated (left) and dry (right) conditions. The bimodal distribution corresponds 
to the occurrence of two different lithologies in the borehole (see also figure [7]). 
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Fig. 2. Calibration of bulk density measurements. The absorption coefficient /i mea- 
sured on dry and saturated samples is converted into bulk density using equation [TJ 
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Fig. 3. Variation of sonic velocity with confining pressure. Top: Laboratory data 
together with Unear-exponential and exponential fit. The quality of the fit differs 
only in the extrapolated range. Bottom: Data before and after pressure correction 
based on equation [3l 
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Fig. 4. Cross plot of properties measured in the laboratory: Thermal conductivity 
A versus sonic velocity Vp (top left), bulk density pb (top right), and porosity cj) 
(bottom right), (p is derived from dry and saturated gamma density measurements. 
Sonic velocity versus bulk density is plotted in the lower left. Black and grey symbols 
correspond to sandy and carbonaceous samples. Dry and saturated measurements 
are shown as (+) and (0), respectively. 
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Fig. 5. Comparison of different mixing laws for thermal conductivity: HS^ - Hash- 
in-Shtrikman upper and lower bounds. HS(avg) - mean of HS"'' and HS" bounds. 
Results are shown for saturated and dry rock samples. Matrix, water, and air ther- 
mal conductivities are taken to be 5W(mK)"^, 0.6W(mK)"^, and 0.026 W(mK)"^, 
respectively. Choosing an inappropriate mixing law results in much larger errors for 
dry samples. 
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Fig. 6. Variation of aspect ratio a versus matrix thermal conductivity Am- Both 
properties are computed using equations [13] and [HI We used dry and saturated 
thermal conductivity and porosity derived from bulk density as input to the model. 
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Fig. 7. Composite log showing caliper, sonic, neutron porosity, gamma-ray, and core 
lithology. 
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Fig. 8. Sand-Shale-Carbonate model for the borehole studied. Input data (left 
panel): NPHI - neutron porosity; GR - natural gamma radiation; DT - acoustic slow- 
ness. Input logs are used to compute the composition (middle left panel). Colour 
coding: black - shale; light grey - sand; dark grey - carbonates; white - porosity. 
Composition is then used to compute saturated and dry thermal conductivity, TCS 
and TCD, respectively (middle right and right panel). Thermal conductivity derived 

from logs is compared to saturated and dry core data, CTCS and CTCD. 
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Fig. 9. Sand-Shale model for the borehole studied. Input data (left panel): GR - 
natural gamma radiation; DT - acoustic slowness. Input logs are used to compute 
the composition (middle left panel). Colour coding: black - shale; light grey - sand; 
white - porosity. Composition is then used to compute saturated and dry thermal 
conductivity, TCS and TCD, respectively (middle right and right panel). Thermal 
conductivity derived from logs is compared to saturated and dry core data, CTCS 
and CTCD. 29 



